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Abstract 

The momentum-dependent local-ansatz (MLA) wavefunction describes well correlated electrons 
in solids in both the weak and strong interaction regimes. In order to apply the theory to the 
realistic system, we have extended the MLA to the first-principles version using the tight-binding 
LDA-i-U Hamiltonian. We demonstrate for the paramagnetic Fe that the first-principles MLA can 
describe a reasonable correlation energy gain and suppression of charge fluctuations due to elec¬ 
tron correlations. Furthermore, we show that the MLA yields a distinct momentum dependence 
of the momentum distribution, and thus improves the Gutzwiller wavefunction. 

Keywords', variational method, local-ansatz wavefunction, Gutzwiller wavefunction, electron cor¬ 
relations 

1 Introduction 

Electron correlations play an important role on the properties of materials such as the magnetism, 
the metal-insulator transition, and the high-temperature superconductivity. The wavefunction 
method is a simple and intuitive method to describe the correlated electrons leading to these 
properties. In this method, we construct a trial wavefunction with variational parameters choosing 
the minimum basis set in the Hilbert space to describe the correlated motion of electrons, and find 
the best wavefunction at the ground state on the basis of the variational principle. 

The Gutzwiller wavefunction (GW) [[LO is one of the oldest and useful trial wavefunctions 
in solids. The wavefunction describes the correlated electrons controlling the amplitudes of the 
doubly occupied states in the Hartree-Fock wavefunction. The GW has been applied to many 
phenomena such as the correlation effects on the magnetism, the heavy-fermion behavior, and 
the metal-insulator transition. The local-ansatz wavefunction (LA) [|l]-0 is an alternative wave- 
function. It takes into account the correlated states expanded in terms of the residual Coulomb 
interactions, which are neglected in the Hartree-Fock wavefunction. The LA has been applied to 
many systems such as the transition metals and the semiconductors [|71[8l . 

The wavefunctions mentioned above, however, do not yield the exact results of the Rayleigh- 
Schrddinger perturbation theory in the weak Coulomb interaction limit. We therefore recently 
proposed the momentum-dependent local ansatz (MLA) wavefunction [l9l-[T^. In the MLA, we 
first expand the Hilbert space from the Hartree-Fock ground state |0o) by means of the two- 
particle excitation operators in the momentum representation. After introducing the momentum 
dependent amplitudes for the excitation operators as the variational parameters, we project these 
operators onto the local orbitals. Resulting correlators {Oi} are more suitable than those created 
by the LA, for the description of electron correlations in solids. We also emphasize that the ex¬ 
act wavefunction in the weak interaction limit is essential to describe the quasiparticle weight in 
the strongly correlated regime, as shown in our paper IITOl . because the latter is obtained by the 
renormalization of that in the weak interaction limit in the variational approach. In the second 
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paper uni, we extended the MLA to the strongly eorrelated regime introdueing a hybrid wave- 
funetion as a starting wavefunetion. The potential of the hybrid wavefunetion ean ehange from 
the Hartree-Foek type to the alloy-analogy type by varying a weighting faetor from zero to one, so 
that we ean start from the best wavefunetion for a given Coulomb interaetion strength. We have 
shown on the basis of numerieal ealeulations that the MLA improves the LA and the GW lllOllllll . 

In this paper, we propose the first-prineiples MLA eombining the MLA with the first-prineiples 
tight-binding Hamiltonian, ealled the LDA (loeal-density approximation)-!-!! Hamiltonian. The 
first-principles MLA allows us to describe the correlations of the multiband system, especially, 
the on-site charge-charge correlations, as well as the Hund-rule correlations at the ground states. 
We limit ourself in the present paper to the MLA in the metallic regime starting from the Hartree- 
Foek wavefunetion. The description of the strongly correlated regime with use of a hybrid wave- 
function will be presented elsewhere. 

In the following section, we propose the first-principles MLA wavefunetion, and obtain the 
correlation energy in the single-site approximation (SSA). Using the variational principle, we 
obtain an approximate form of the variational parameters. In Sec. 3, we present some numerical 
results of calculations for paramagnetic Fe in the weak Coulomb interaction regime. We verify the 
correlation energy gain, the suppression of charge fluctuations, and in particular the momentum 
dependent behavior of the momentum distribution. In the last section 4, we summarize the present 
work. 

2 LDA-fU Hamiltonian and MLA Wavefunetion 

We consider here the transition metal system for simplicity and adopt the first-principles LDA-i-U 
Hamiltonian, which is based on the tight-binding linear muffin-tin orbital method [[T3l[T4ll . 

H = Ho + H,, (1) 

^0 = ALa + XI (2) 

iLcr iLjL'a 




L Umm H" 
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' m>m' 


(3) 


Here Hq and Hi denote the non-interacting and interacting parts of the Hamiltonian H, re¬ 
spectively. is the atomic level of the orbital L on site i. tnjL' is the transfer integral 
between iL and jL'. L = {l,m) denotes the s,p,d orbitals. al];^^{aiLa) is the creation (an¬ 
nihilation) operator for an electron on site i with orbital L and spin a, and nua = alj^^^ana 
i^iL = X]a 7 ^lLai^)a'y IS the chaTgc (spin) density operator. Here cr denote the Pauli spin 

matrices. 

In the LDA-i-U Hamiltonian we assume that the sp electrons are well described by the LDA 
in the band theory, and take into account only on-site Coulomb interactions between d{l = 2) 
electrons in the interacting part Hi. Umm{Umm') and Jmm' denote the intra-orbital (inter-orbital) 
Coulomb and exchange interactions, respectively, humisum) with I = 2 is the charge (spin) 
density operator for d electrons on site i and orbital m. The atomic level in Hq is calculated 
from the LDA atomic level 6*^ by subtracting the double counting potential as — 

9Ei^iQjs^/dniLa- Here UiLa is the charge density at the ground state, i^LOA L a LDA functional for 
the intra-atomic Coulomb interactions. 

In the first-principles MLA, we rewrite the Hamiltonian H as the sum of the Hartree-Foek 
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Hamiltonian Hup and the residual interactions Hi. The latter is given by 
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Here 0-mm’ ^imm” ^imm' denote the two-particle operators: the intra-orbital operators, the 
charge-charge inter-orbital operators, and the spin-spin inter-orbital operators, respectively, and 
are defined as follows. 






(5) 
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Here Snumcr ’^Uma (j^ilmcr)Ot ^ilm {'^ilm}oj and ^ilm (^iZm)oj (~)o 

denotes the average in the Hartree-Fock approximation. 

In the residual interactions dH), we have three types of correlation operators, 

( 2 ) 

and Accordingly, we construct the MLA wavefunction as follows. 
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The correlators O-mm’ ^^d *^he two-particle operators projected onto the local 

orbitals, and describe the intra-orbital correlations, the charge-charge inter-orbital correlations, 
and the spin-spin inter-orbital correlations, respectively. They are defined by 


O'iL' = (k 2 n 2 liL}(iLjk 2 n 2 )(k[n[jiL'}(iL'jkini) 

{kna} 


(9) 


Here q =0, 1, and 2. al.^^{akna) is the creation (annihilation) operator for an electron with 
momentum k, band index n, and spin a. They are given by those in the site representation as 

^kna r 0^iLa-{k'kl\iL). 

The momentum dependent amplitudes ^pl/^ 2 ' 2 i’i} given by 

'^LL'{2'21'1} VLk'2n'2k2n2k[n'j^kini^LL'Sa!2i^a-2"[^o-[l^cri'[) 
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Here 'Z7Lfc2’*2^2fi2fc(n(fcini5 Cii'fc'fi,^fc2n2fc(n(fcini’ '^LL'fc2^l■2^2^^2Zc(n(A:lr^l ’ ^^d 'b2y2y'fc2^^2^2^^2Zc(n(fclr^l 
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It should be noted that the MLA wavefunction ([8]) describes exactly the weak interaction 
regime. The best wavefunction is obtained by choosing the variational parameters best on the 
basis of the variational principles for the ground state energy. 

The correlation energy is defined by (if) — (if)o, where ( )o) denotes the full (Hartree- 

Fock) average. The MLA contains an infinite number of variational parameters, so that the ad¬ 
vanced numerical methods such as the variational Monte-Carlo technique are not applicable for 
the calculation of the energy. We adopt here the single-site approximation (SSA). The correlation 
energy per atom ec is then obtained as follows. 


-(o/^)o-(^0,)o + (o/^0,)o 
1 -f {Pi Oj)o 


Here H = H — (if )o, and 
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Solving the self-consistent equations obtained from the stationary condition 5ec = 0, we 
find an approximate form of variational parameters in the weak Coulomb interaction regime as 
follows. 
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Here fciniCTi *^fc 2 fi 20'2 TCfc^njo-'i ^k\n\cr\ dcnotcs the two particle ex 

citation energy. i)m, are the variational parameters. Substituting Eqs.(fT3l)~(fT^ 

into Eq. lfTT]) and using the variational principles 5ec = 0 again, we obtain the self-consistent 


equations for the variational parameters to be solved. The elements {Oi ff)o, {Oi HOi)^, and 


(o/di)o in Eq.lfTTl) are calculated by Wick’s theorem. 

The other physical quantities such as the electron number and charge fluctuations are also 
obtained by taking average with respect to the MEA wavefunction ([8]) and making the SSA. 


3 Numerical Results in the Weak Interaction Regime 

The ferromagnetism and related properties of Ee have extensively been investigated theoretically 
with use of realistic Hamiltonians with s, p, and d orbitals at the ground states IITSlfT^ and 
at finite temperatures IfTTUT^ . But its physics has not yet been fully clarified. We performed 
numerical calculations for the paramagnetic Ee in order to clarify the basic behavior of the first- 
principles MEA in the weak Coulomb interaction regime. In the weak Coulomb and exchange 
energy limit, we obtain fjm ^ I, Cmm'^ 1’ ^mi' -1, and -)■ -1 from the self- 

consistent equations for the variational parameters. We adopted in the present calculations these 
variational parameters, and calculated physical quantities of Ee. Eurthermore we neglected the 
orbital dependence in the Coulomb and exchange integrals in the numerical calculations, and 
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adopted the values Umm = Uq = 0.2749 Ry, Umm' = Ui = 0.1426 Ry, and Jmm' = J = 0.0662 
Ry which are obtained from the relations Uq = U + 8J/5,Ui = 17 — 2J/5, and J = J, using the 
average values U = 0.1691 Ry and J = 0.0662 Ry by Anisimov et al. [|20l . Note that we adopted 
here the relation Uq = Ui + 2J for cubic system. The transfer integrals and the atomic level have 
been calculated with use of the Stuttgart tight-binding LMTO (linear muffin-tin orbital) package. 



Figure 1: The correlation energy Cc as a function of Coulomb interaction strength alio for the 
paramagnetic Fe. Dashed curve: the result without inter-orbital correlations (i.e.,Ui — Jl2 = 
J = 0), solid curve: the result with both the intra- and inter-orbital correlations. 



Figure 2: The charge fluctuation {{SnidY) vs Coulomb interaction strength aUo curve for the 
paramagnetic Fe. Dashed curve: the result without inter-orbital correlations, solid curve: the 
result with both the intra- and inter-orbital correlations. 


In order to see a systematic change due to correlation strength, we scaled Uq, Ui, and J as 
aCo, aUi, and aJ using a scaling factor a. Figure [H shows the calculated correlation energy as 
a function of aUo. With increasing aUo (as well as aUi and aJ), we find that the correlation 
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energy ec monotonieally deerease and find ec = —0.0516 Ry when a = 1 and Ui = J = 0. When 
we take into aeeount the inter-site eorrelations, the eorrelation energy €c deereases further and we 
obtain €c = —0.1101 Ry when a = 1. 

The eorrelation energy gain is aeeompanied by the suppression of eharge fluetuations. We 
ealeulated the eharge fluctuations for d electrons {{SuidY) = (nj^) — {nidY ^ function of alio. 
As shown in Fig.[2l the charge fluctuation in the Hartree-Fock approximation is 2.2. The intra¬ 
orbital correlations suppress the charge fluctuations and yields ((5nid)^)=1.73 for alio = 0.2 Ry. 
The inter-orbital correlations more rapidly decrease the charge fluctuation with increasing alio as 
seen in Fig.[2l Calculated charge fluctuation is {{SriidY) ~ 1-3 for alio = 0.2 Ry. The result is 
comparable to the value of the LA with the d-band model |I71, {{SriidY) ~ 1-0, but is somewhat 
larger than that of the LA because the present theory takes into account the hybridization between 
the d and sp electrons. 
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Figure 3: The momentum distribution {riknY as a function of the energy ikna for various scaling 
factors a of the Coulomb and exchange energies. Dash line denotes the distribution in the Hartree- 
Fock approximation. 



We emphasized in the previous paper for the single-band model lIT^ that the MLA can de¬ 
scribe the momentum dependence of the momentum distribution, while those calculated from the 
original LA and the GW remain almost constant below and above the Fermi level. This behavior 
is seen also in the first-principles case. In fact, the momentum distribution {nkna) is given in the 
present theory as follows. 


/ \ _ / \ I N{Qi nknaOi)o 

Y^kna) {^kna/O _ 4 - _ • ( 1 ^) 

1 + {Oi Oi)o 

The first term at the rhs (right-hand-side) is the momentum distribution in the Hartree-Fock ap¬ 
proximation, which is given by the Fermi distribution function at zero temperature /{ekna) = 
d{—ekna)- Here 9 denotes the step function, and ikna is the Hartree-Fock one-electron energy 
measured from the Fermi level. The second term at the rhs of Eq. (fTTI) describes the correlation 
correction. N denotes the number of atoms, hkna is defined by fikna = nkna — {nknYo- The 
correlation correction consists of the terms being proportional to \uLna-{k)\'^ fYkna) and those 
being proportional to \uLnaik)\'^ fi—^kna)- Here {umaik)} are the eigenvectors for given k 
points. If we replace |MLncr(fc)P with 1/5 as a rough approximation for the d-like branch n near 
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the Fermi level, we find that (nkna) depends on the momentum only via the energy ikna as in 
the single-band model. Figure [3] shows the calculated result of momentum distributions in this 
approximation. We find clearly the momentum dependence of (rikno) via iknu, and obtain the 
mass enhancement m* = 1.4 for a = 1.0 from a jump at the Fermi level. This value should 
be compared with the experimental renormalization values 1.38 ~ 2.12 which are obtained from 
the comparison of the LDA band calculations with the T-linear specific heat at low tempera¬ 
tures [|15[I21II . We have to calculate {ukna) along the high-symmetry lines taking into account the 
fc-dependence of luLnaik)]"^ in more detailed calculations. The results of the full calculations for 
{f^kna) will be published elsewhere. 

4 Summary 

We have extended the MLA to the first-principles MLA on the basis of the tight-binding LDA 
-i-U Hamiltonian. The wavefunction is constructed by applying the momentum-dependent intra¬ 
orbital correlators, the inter-orbital charge-charge correlators, and the spin-spin correlators to 
the Hartree-Fock uncorrelated wavefunction. The theory yields the exact results in the weak 
Coulomb interaction limit, and describes the charge-charge correlations as well as the Hund rule 
correlations in the real system. 

By means of the numerical calculations for the paramagnetic Fe in the weak Coulomb in¬ 
teraction regime, we verified the correlation energy gain due to the inter-orbital charge-charge 
correlations as well as the spin-spin correlations between the d orbitals. The charge fluctuation 
{{dnidY) in the present calculations is 1.3. The value is somewhat larger than the value obtained 
by the LA with the d-band model because of the hybridization between sp and d orbitals. We 
also clarified that the present theory leads to a distinct momentum dependence of the momen¬ 
tum distribution. The result qualitatively differs from those obtained from the LA and GW. The 
calculated effective mass m* = 1.4 is consistent with the experimental values obtained from the 
low-temperature specific heat data. Numerical calculations in the present paper are limited to the 
weak Coulomb interaction regime. In order to examine the quantitative aspects of the properties 
of correlated electrons in Fe and its compounds, we have to solve the full self-consistent equa¬ 
tions for the variational parameters {i.e., pm, Numerical calculations for 

the correlated electron regime are in progress. 
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